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The waiting time distribution w(t), i.e. the probability for a delay r between two subsequent 
transition ('jumps') of particles, is a statistical tool in (quantum) transport. Using generalized 
Master equations for systems coupled to external particle reservoirs, one can establish relations 
between w(t) and other statistical transport quantities such as the noise spectrum and the Full 
Counting Statistics. It turns out that w(t) usually contains additional information on system 
parameters and properties such as quantum coherence, the number of internal states, or the entropy 
of the current channels that participate in transport. 
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I. INTRODUCTION 

Particle transport through a given open system S can 
be regarded as a form of spectroscopy: by monitoring the 
flow of particles between the system and external particle 
reservoirs (usually in a stationary state), one likes to ex- 
tract as much information about S as possible. Take S as 
a 'black box' and the series of times t"(i — 1, 2, ...), where 
individual particles leave or enter S through reservoir a, 
as the only available data. Just by counting (and per- 
haps labeling according to energy, spin etc.), one would 
like to 'reconstruct' S, thus defining an inverse problem 
where (some) transport data are known but the system 
Hamiltonian or Liouvillian is not. 

Of course, part of this scenario is just a description of 
modern (quantum) transport, i.e. the determination of 
full counting statistics (FCS) or higher cumulants, cur- 
rent fluctuations, or noise spectra S(u>), rather than the 
determination of, say, 'just' the stationary current or a 
current- voltage characteristics. 

In this paper, the waiting time distribution w(t) as 
a statistical tool for transport (described by generalized 
Master equations) is analyzed for systems with typically 
a finite number of discrete internal states (including pos- 
sible quantum coherences), among which transitions oc- 
cur due to the coupling to external reservoirs. w(t) is 
a probability density for the delay time r between two 
subsequent 'jump' events (e.g. single electrons tunneling 
out of a quantum dot) and is well-known from quantum 
optics^. In quantum transport, it has occasionally been 
used in the discussion of shot noised, counting statistics^ 
and in the description of transport through single, vibrat- 
ing molecules^. The purpose of this paper is to find out 
in how far the waiting time distribution contains informa- 
tion that is (partly) complementary to the one contained 
in, e.g., FCS and S(ui). For example, these latter can be 
derived from w(t) but not vice- versa. 

In particular, the waiting times contain spatially sep- 
arate information on, e.g., tunnel barriers in quantum 
dots, and they depend on the number of internal transi- 
tions, the quantum coherence, and the single or multiple 
'reset' character of the system, i.e. the number of system 
states that are directly involved in the process of parti- 



cles entering or leaving the system. The simplest example 
of a single-reset system in quantum optics is resonance 
fluorescence in two-level atoms, where each emitted pho- 
ton leaves the atom in its ground state. Corresponding 
examples in electron transport are quantum dots in the 
regime of strong Coulomb blockade, where only tunnel- 
ing of one additional 'transport' electron is possible. As 
a consequence, w(t) then has to vanish for r — > when 
measured at the same terminal. 

Experimentally, various groups-»2> have achieved to 
measure time-series of single electron tunneling events, 
e.g. by monitoring transport through quantum dots us- 
ing a nearby quantum point contact. One has to point 
out, however, that the formalism developed here always 
refers to counting the additional (or missing) particles in 
the reservoirs and not in the system S. Counting transi- 
tions within S by direct interaction of a counting (mea- 
suring device) with S usually would destroy quantum 
coherent features within S. 

The paper is organized as follows: the next section 
introduces the general method on a more formal level 
which is introduced for generalized Master equations. In 
the third section, various example applications and the 
entropy of currents as a possible way to distinguish be- 
tween different classes of systems are discussed. Some 
remarks are left for the conclusions. 



II. METHOD 

There is a large number of cases where one can describe 
quantum transport in terms of a Markovian generalized 
Master equation for the reduced system density operator 
P(t), 

M 

P (t) = Cp(t), £ = £ + A, £i = J2 J >f C 1 ) 

k=l 

Here, L\ describes M different types of jump processes 
where single (quasi) particles tunnel into or out of the 
system. Most of the applications below refer to single 
electron tunneling, but the formalism stays valid as long 
as one has a time evolution of the form Eq. fTJ). 
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The individual jump processes depend on the specific 
terminal (e.g., left and right), energy level, system state 
before and after the jump etc., but at this stage one 
makes no further assumption other than that the jump 
operators act on density matrices p as 



JkP = c k [p]pk, 



(2) 



where pk is a density matrix and Ck[p] is a scalar that is 
typically given in terms of positive transition rates. 

The splitting C — Co + C\ is not unique and depends 
very much on the kind of information one wishes to ob- 
tain from a transport measurement. A large M means 
that one can access a large number of jump processes of 
different kind (e.g., energy resolved), but this is typically 
very difficult to monitor experimentally. Also numeri- 
cally, there is a trade-off between making Co as diagonal 
as possible by having a large M, leading to large but 
simply structured matrices for the waiting times and a 
corresponding large inversion problem for the noise spec- 
trum, and having small M, leading to small waiting time 
matrices with, however, a relatively complicated struc- 
ture (see below). 

In the interaction picture with respect to Co and 



St 



o C t 



a formal solution ('unraveling') of the Mas- 



ter equation Eq. ([T]) reads 



= e e / / ' 

n=Oh=l,...,l n =l Jo Jo 



(3) 



where the 



PL 



Ji ({**}) = St-t n Jl n St n -t n - 1 Jl n _ 1 ---Jl 1 St 1 Pin (4) 



are (unnormalized) conditioned density operators^ that 
describe the non-unitary time-evolution of an initial den- 
sity operator p ln , interrupted by n quantum jumps of 
type U at times U with i = 1, ...,n. 

In the following, the existence of a unique stationary 
solution po of Eq. (TT|) is always assumed, i.e. decompos- 
able Liouvillians^ with a block structure referring to the 
uninteresting case of two (or more) de-coupled systems 
are excluded. For n = 2 subsequent jumps of type k and 
I with t\ = and t% = r = t, one now defines the nor- 
malized conditioned density operator pfa with p ln = po 
inEq. gD, 



P C ki(r) 



J k e CoT Jip 
w k i(T)TiJip a ' 



(5) 



The normalization factor in the denominator defines a 
matrix W(r) of waiting time distributions 



w u (r) = (W(r)) 



ki 



where 



h = TrJipo 



(6) 



(7) 



defines a stationary current due to jump processes of type 
I. Here and in the following, all currents have physical 
dimension 1/time. 

Eq. fl} defines a stochastic process that is described 
by a linear system of coupled first order differen- 
tial equations. For practical calculations, it is there- 
fore convenient to represent the density operator p 
as a column vector with real entries, for example as 
p= {pn,p 2 2,--,PNN,'3pi2,'3tpi2,--,'$tpNN-i) T - Corre- 
spondingly, the super-operators J7fc, C, Co become real 
matrices. The next step is to introduce a convenient 
Dirac-like notation^ where kets p <-> \p)) denote nor- 
malized density operators. The stationary state po with 
Cpo = is denoted as ket |0)). The ((0| is the row vector 
(1, 1, ...1, 0, 0...0) such that in this vector representation 
of density operators p, the trace operation on a column 
vector Ap (where A is an arbitrary super-operator) be- 
comes the scalar product TiAp = ((0|^4p)). In particular, 
one has ((0|p)) = 1 for density operators p. Finally, the 
bras ((k\ are real row- vectors defined via the action of the 
jump operators J7fc, Eq. when writing their matrices 
as the dyadic product 



Jk = 



(8) 



Note that bras (row vectors) and kets (column vec- 
tors) |fc)) are independent vectors and are not dual to 
each other. They are introduced here as a convenient no- 
tation for the calculations to follow. They have positive 
entries as the transitions rates for the quantum jumps 
are positive. 

As the diagonal elements of the Liouvillian C do not 
describe jump processes but the conservation of proba- 
bility, the matrices J k have no diagonal elements, and 
consequently ((k\k)) — for k ^ and J% = 0. 
By furthermore taking the trace in Eq. @, one has 
c k [p] = TiJ k p = ((k\p)} and also in Eq. ©, 



Wkl (t) 



((k\e^\l))((l\0)) 

«i|o» 



«fc|e^|Z», (9) 



where the stationary currents // = ((^|0)) due to jump 
processes of type I cancel, cf. the definition Eq. ©. As a 
consequence, the calculation of the waiting times wu (r) 
does not require the knowledge of the stationary state, po- 
Note that the wh{t) are manifestly real quantities. With 
((k\ and having positive entries, the Wki{r) are also 
positive, provided that the time-evolution operation e c ° T 
does not lead to non-positive definite density operators 
- a condition which is physically plausible and which is 
fulfilled in all the examples discussed below. 

One now defines a waiting time super-operator Wi(z) 
as the Laplace transform of e c ° T Ji, 



Wi{z) = {z-C )- l J h 



(10) 



by which the Laplace transform Wki(z) of the waiting 
time distributions can be conveniently written as an 'ex- 
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pectation value'. 

wm(z) = 



dte- zt w k i{t) 



TvJ k Wi(z)p 



= miz-car'ii)). 



(ii) 



For any initial jump of type /, the waiting time dis- 
tribution Wki (i) has to give unity when summed over all 
subsequent jumps of type k and integrated over all times 
t. This normalization indeed follows from 



/>OC 

/ dtJ2w k i(t) =J2w kl (o) = 

J ° k k 

((bKc-co^jtio)) = mm) 



(12) 



1, 



where one exploits the fact that ((6| is a left eigenvector 
with eigenvalue zero of the total Liouvillian, ((0|£ = 0, 
which reflects conservation of probability. 



A. Relation to the Noise Spectrum 

(Quantum) fluctuations of the system can be visualised 
in fluctuations of the current and play a central role 
in analysing the internal system dynamics via transport 



spectroscopy 



12.13.14.15 



S k i(u) = - 



The spectrum 
dte iw *({57 fc (t),^(0)}) 



(13) 



defines fluctuations of currents, where 8Ik(t) denotes the 
deviation of current Ik from its average in the station- 
ary state, Eq. {?]). It is derived via the MacDonald 
formula 1 ^ 16 !^, 



Ski(v) 



Sail 



(14) 



2 I ±iu + £ 



+ (fc<->0 



which for numerical convenience can also be expressed 
in terms of a resolvent operator 11 R(uj) = i\0))((0\/u + 
Q[iuj + C}- 1 Q with Q = 1 - |0))((6| and where the sin- 
gular contribution from the stationary solution at lo = 
has been projected out. 

The link with the waiting time distributions is now 
established via the operator identity 



(z-C)- 1 ^ = \{z-C )(l-{z-C )- 1 d 



A I 



W(z). 



ji 

(15) 



The expression Jk (z — C) Ji is then formally ex- 
panded into a geometric series, the first term of which 
is given by 



({0\J k Wi(z)\0}) 



iB w (*)/j = /i[W(*)] w , (16) 



where W(z) is the matrix of the Laplace transformed 
Wki (z). The n + 1-th term of the geometric series contains 
correspondingly 

((k\W mi (z)...W mn {z)Wi(z)\0)) = (17) 

= W kmi (z)w mim2 (z)...W mn l(z)Il (18) 

which upon summation simply yields matrix products, 
i.e. the n + 1-th power of the matrix W, and thus 



{{Q\J k {z- Cy 1 Jt\Q)) 



(19) 



M 



kl 



= (WE E w »( z > mm) 

n=0 \m=l / 

oo 

= J2 h [W" +1 (z)] kl = Ii \(1 W(z))' 1 W(z) 



Up to here, this is a formal result as the sum Eq. (Tit 
does not need to converge (the expression actually has 
to be singular for z = due to the normalization of the 
waiting times). For purely imaginary z = iu> with real 
lu =/= 0, however, for positive Wki(t) > and due to the 
normalization Eq. (|12p . the matrix norm 

||W(w)|| = max^ \ w k i(ioj)\ 

k 

poo 

< max V / dtw k i(t) = 1, (20) 
1 k Jo 

and with ||W"(z)|| < ||W(z)|| n it follows that the series 
converges. 

Purely imaginary z — iuj is just what is needed in Eq. 
(|T4"|). and one obtains 



Ski{u) = hill 



(21) 



\ I E [c 1 - w^r 1 w(±»w)] u h + (k <-> i) 1 , 



which expresses the noise spectrum in terms of the wait- 
ing time distributions. 



B. Relaxation Currents 

After a quantum jump of type I, the system relaxes 
from the reset state pi = \l)) into the stationary state. 
Such a relaxation potentially involves all jump processes 
of type k, and consequently relaxation currents should 
be defined as 

r k f™(t) ee TrJ k e ct Pl = «fc|e £t |0). (22) 
Upon Laplace-transformation, one finds 

TYjk(z-c)- 1 j lPo = ((^(z-rrMOXailo)) 

fre 
'kl 



Ilf ax (z)I h z^O, (23) 
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where P h f ax (z) is the Laplace transform of (^ lax )(i), and 
thus by comparison with Eq. (114[) , 



SM = -iA5 kl 



l (±iuj) } + (k 0( 24 ) 



The fluctuations of the stationary currents are thus de- 
termined by the relaxation currents from the reset states 
pi. The relaxation currents are in turn determined by 
the waiting time distribution, 



r k f^(z)= (i-w(«))- x w(2) 



(25) 



In the time domain, the relaxation currents are obtained 
from the waiting times by re-expanding Eq. (|25p as a 
series of convolution integrals, 



w k i (t) 



f «fti[W(t-*i)W(ii)] 
Jo 



H 



(26) 



Eqs. (|2"5)) also illustrate a numerical trade-off 

when calculating waiting time distributions and noise 
spectra for different unravelings of £, Eq. ([T]): splitting a 
large number M of jump operators off the total Liouvil- 
lian C has the benefit of potentially simple Liouvillians 
Cq, e.g. close to lower-triangular form. As a downside, 
this leads to large waiting time matrices W(z) and there- 
fore to a potentially large matrix inversion problem for 
the determination of the If.f &K (z) when calculating the 
noise spectrum. 



C. Single and Multiple Reset Systems 

The individual jump operators Jk, Eq. ([2]), refer to a 
unique density matrix |fc)) after the jump of type k. One 
can now further classify jump processes by introducing 
class labels a by slightly extending the notation, writing 
double indices ka in 



£ — Co + £i, 



£i — ^ Jka 

ka. 



(27) 



with Jk a = \ka))((ka\, cf. Eq. @. In the examples 
below, a = L/R labels two (left and right) particle reser- 
voirs (leads) attached to the system. The index k then 
further specifies the type of jump process for reservoir 
a. Depending on the physical system described by C, k 
can label the jump process according to energy, spin, or 
some other degrees of freedom that are coupled to the 
transport process, e.g. phonons. 

A simple counting without energy resolution of, e.g., 
the emitted electrons, is the situation which - in contrast 
to quantum optics - is typical of quantum transport ex- 
periments. One then is interested in sums of jump op- 
erators Jka only, i.e. super-operators J a and currents 

la 



E 

k 



Jka, la = ^JaPO- 



(28) 



Individual jump processes within the class a typically 
leave the system in different states \ka)). Some systems 
have a unique 'reset' state |a)), as for example the empty 
state of a quantum dot in the strong Coulomb blockade 
regime with only one or zero electrons. This is in con- 
trast to other cases where for example the system parti- 
cle number fluctuates by more than one in the stationary 
state. In general, this distinction leads to the definition 
of single and multiple 'reset' systems: one defines a sin- 
gle reset system as a system where all jump operators 
J a = Ylk ^ ka nave a reset property analogous to Eq. (J2J), 
i.e. 



JaP = C a [p]p a . 



(29) 



This condition is fulfilled if the state \ka)) = \a)) inde- 
pendent of k, i.e. 

J a = J2 \ka)){{ka\ = |a» ^{{ka\ ee \a))((a\. (30) 



In general, a sum over jump operators Jk can not be 
written in 'separable' form Eq. (f3"0")) . but the bras and kets 
are 'entangled', and the corresponding system is called a 
'multiple reset' system. 

In analogy with the definition Eq. (fTTj) of the waiting 
time distribution, the quantity 



TrJ a Wp(z)p Q 
TrJpPo 



(31) 



defines a (reduced) waiting time distribution which in- 
volves the J a , Eq. (|2"8")l . only. Accordingly, one defines 
the reduced noise spectrum matrix of the currents I a and 

Iff, 

S^(u) = S aP I p (32) 

= \h {<w + E ^t x (±^)| + 

with S^l(uj) = J2ki Ska.ipi^) an d the reduced relaxation 
currents 



frclax 
1 a/3 



Em [(l-woor'wc* 



h 



ka,l/3 



Eiiw 



.(33) 



Analysing this definition, one recognizes a profound dif- 
ference between single and multiple reset systems in their 
respective connections between noise and waiting time 
distribution: in single reset systems, due to the unique- 
ness \lfi)) = one has 



{{a\{z-C a y l \P)). 



(34) 
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The reduced relaxation currents then simplify drastically: 
re-expanding the geometric series in Eq. (|33[) . 

k 

+ EE«H( Z - £o)- V»((^ - A.)" 1 !/?)) + ... 



k I'. 



where the sum over Z cancels the currents X^-fyS) 
therefore 



(35) 



and 



'a/3 



(l-W (r) (z)) W (r) (^) 



. (36) 



a/3 



with the matrix W^'^(z) of the reduced waiting time dis- 
tribution, Eq. (f3"Tj) . Therefore in single reset systems, 
the reduced noise spectrum matrix S^I(u>) and the re- 
duced waiting time distribution w^l(ui) are connected 
via a matrix relation, cf. Eq. (|36p and Eq. (|32p . involv- 
ing the matrix ~W^(z) which in general has a smaller 
size than the full waiting time matrix W(z). 

In the non-separable case, i.e. where the super- 
operator for the total current through the system can 
not be written as in Eq. ([30]) . this is no longer the case, 
and one has to use the full matrix equation Eq. (f2"Tj) in 
order to relate both quantities. 



One obtains the short time expansion of wu (t) from an 
expansion of Wki(z) into a Laurent series. For example, 
when 



Wki{r) ~ c k iT 71 







(40) 



with n > 0, one has Wki(z — > oo) ~ CkiT(n)z~ n , where 
r(.) is the Gamma function. Alternatively, the exponent 
n is directly obtained by using the expansion 



{{k\e^\l)) = 

{{k\l + £ a T+ T -Cl + ...\l)). (41) 



In the examples below, one finds that the asymptotic be- 
haviour of Wki(r) at small r, i.e. the exponent n, depends 
on the dynamics within the system and thus contains po- 
tentially useful information. 

Finally, the high frequency expansion of the reduced 
noise spectrum, Eq. (|3"2|) , reads 




(±iuj) 



(42) 



which means that at large frequencies, the correction to 
Sj1(o)) to the Poissonian limit is directly given by the 
reduced waiting time distribution regardless of whether 
or not the system is of single reset type. 



D. High-frequency and short time expansions 



E. Relation to Full Counting Statistics (FCS) 



The noise formula, Eq. (f2"Tj) . can formally be expanded 
in powers of Wki(z), 



iL0) + O(w 2 ) \h 



+ (fc «-►!)> 



(37) 



where the terms involving first and higher powers of 
Wki(z) describe the deviation from the Poissonian limit 
Sm(lo) = Skill- In actual fact, this limit is reached for 
lu — > oo since liniz—yoo Wki(z) = 0, and one therefore ex- 
pects Eq. (|37p to be a high-frequency expansion. In par- 
ticular, one has for z — > oo, 



ilf^(z)-,w kl (z), z 



(38) 



The large z limit corresponds to small waiting times r, 
and for small r the system has not enough time to re- 
solve all its possible states after a quantum jump. At 
small times i, relaxation currents and waiting times thus 
coincide, 



JI.f ax (r) -> w m (t), t 



0. 



(39) 



which also follows from Eq. 



For any given unraveling C = Cq + C%, Eq. (JTJ) , one 
expects the waiting time distribution to be related to 
the Full Counting Statistics (FCS), i.e. the probability 
p(n,t) of n quantum jumps in a time interval [0,t]. For 
a single reset system and a single jump operator C\ = 
C = the generating function G(x,t) associated 

with p(n, t) is obtained from the counting-variable (x)- 
dependent propagato r 19 i 20 i 21 i 22 , 

oo 

G( X , t) = J2 e mx p(n, t) = Tre^+^Vo, (43) 

n=0 

where po again denotes the stationary state. The Laplace 
transform G{\, z) of the generating function therefore is 

G( X ,z) = ({OKz-Co-jxj)- 1 ^)) 

= E«6|[(^-A)- 1 e ix J] n (^-£o)- 1 |0)) 

w i(z)w w (z) 



n=0 

woo(z) 



(44) 



with the definition 

w tJ (z) = ((iKz-^o)" 1 !^), i,i,= 0,l. (45) 
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The long-time behaviour of p(n, t) is then obtained from 
the x-dependent pole of Eq. (|4"4"|) . i.e. the solution z of 



(46) 



operators Jk that are included in C\ = Ylk=i ^ k defines 
the jump processes that one wishes to monitor. 

Consider two splittings involving only one and two 
jump operators, 



e lx - w n (z ) = 



with zo(x — 0) = (note that the normalization of t»u 
is Wn(0) = 1), and derivatives of the function Zo(x) then 
yield all the cumulants of p(n, t — > oo) 21 i 22 . 

In the multiple reset case or in general for more than 
one jump operator, one has C\ = YlkLi >^fc> ^1- 0' an d 
correspondingly M counting variables Xk, but the deriva- 
tion of G({xk},z) is analogous to the above one. The 



X-dependent propagator is £({xk}) = A) + YlkLi ^ 
and the generating function is 



G({ X k},z) = w 00 (z) + u T [e-^-W(z)] 1 v, (47) 

where e~ %x is the diagonal matrix of the e~ %Xk and the 
vectors u and v have components 

[ U ] fc = mz-cr^k)) 

[v] k = {{k\{z-C )- l \Q)), k = l,...,M. (48) 
The condition 



det [e~ ix - W(z)] = 



(49) 



then defines a polynomial in z, of which the zero zo({Xfc}) 
with zq(0) determines the FCS. 



III. INFORMATION CONTAINED IN THE 
WAITING TIME DISTRIBUTION, EXAMPLES 



£ = £o + Ji+J 2 =£ + Ji- 



(50) 



The corresponding waiting time distributions are denoted 
as Wij(r) and w(t). Using Ji = the relation 

between these is found via 

= ((lKz-A)- 1 !!)) 

= «li [1- (^-^ )- 1 l2»«2j]- 1 (^-^o)- 1 !!)) 

Wl2(z)W2l{z) 



wn(z) 



1 - w 22 (z) 



(51) 



Two subsequent type-1 jumps are separated by zero, one, 
two,... or any number of type-2 jumps. Therefore, the 
waiting time distribution w(t) for type-1 jumps has to be 
given by the sum over all intermediate jump processes, 
which by re-expanding the geometric series is Eq. (|5ip in 
the time domain, 



w(t) = ion (t) +/ dtiwi 2 (t - ti)w 2 i(ti) (52) 

+ / dti / dt 2 w 12 (t ~ti)w 22 (ti-t 2 )w 21 (t 2 ) + ... 
Jo Jo 



Relations similar to Eq. (|5ip can be derived when more 
than two jump operators are involved. 



From the considerations so far, two questions arise: 
first, what additional information (as compared with the 
noise spectrum Ski(u)) does the waiting time distribution 
wm(t) contain at all. Clearly, the noise spectrum is an 
even function of us whereas w(us) in general is not. In the 
examples below it is shown that indeed w(ui) can reveal 
information that S(u>) does not contain. 

Second, the distinction between various types of jump 
processes gave rise to a matrix structure of both, the noise 
spectrum and the waiting time distribution. When only 
certain combinations of jump operators are accessible ex- 
perimentally (e.g., when electron tunnel is not resolved as 
a function of energy), the difference between single and 
multiple reset systems becomes important and noise and 
waiting time can contain complementary information. 

Before discussing specific examples, two further useful 
concepts are introduced - first, a rule that relates various 
waiting times for different splittings of £, Eq. |T]), and 
second, the entropy of a stationary current distribution. 

A. Waiting Times for Two Different Splittings 

As mentioned above, the splitting of the Liouvillian 
£ = £q + Ci, Eq. ([T]), is not unique. The set of jump 



B. Entropy of Current Distribution 

The currents I a usually split into individual contribu- 
tions Ika = TrJkaPOi for example according to energy or 
spin. A useful concept to quantify the distribution of the 
Ika is the entropy E a of the individual currents Ika, 

E a ^-Y, Ij r^ Ij r- ( 53 ) 

K 

This definition is quite in analogy with the usual en- 
tropy S[pk] = — Pk l°gPfc °f a discrete distribution 
Pk- These entropies, however, depend on the choice of 
the Ik a in the original unraveling. 

In a multiple reset system, a current I a consists of at 
least two individual currents Ika, corresponding to the 
(at least two) reset states \ka)) whence the current en- 
tropy of a multiple reset system is non-zero. In con- 
trast, in single reset systems there is always a choice 
I a = \a))((a\, cf. Eq. (|3H|) . with a single \a)) for which 
the entropies E a are zero. 

In the following, the above concepts are discussed for 
various examples. 
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C. Transitions in a Ring 



contrast, the corresponding noise spectrum, 



Consider a ring with transitions between N + 1 states, 
q 1 2 ... ->iV^0at rates T,, i = 0, 1, ...,N. 
The corresponding Liouvillian is 



C 



/-To ... 

r -ri o 
o r a -r 2 ... 

V o o 



T N \ 





(54) 



For A^ = 1 this is, for example, equivalent to transport of 
single electrons through a single level quantum dot (see 
below) . 

Introducing the single jump operator J = with 
ee (0, ...,T N ) and EE (1,0, ...,0) T , and a splitting 
C = Cq + J corresponding to measuring the jumps N — > 
only, by inverting z — Co one finds the waiting time 
distribution 



){z) = {{l\^-^r 1 \l))=T N [{z~Co)- 1 ] m 
To Ti Tjy 



z + r z + T\ z + r 



JY 



(55) 



where the matrix element of the inverse matrix is a ratio 
of two determinants. Therefore, the waiting time dis- 
tribution is a simple product of N + 1 terms, each of 
which corresponds to a single sequential transition along 
the ring. In the time domain, one has a corresponding 
convolution of exponentials 



U(t) = T ie 



-Tit 



(56) 



which each on its own describe the elementary process of 
independent random transitions. The fact that w(z — > 
oo) ~ To---Tnz~( n+1 ' > leads to the short-time expansion 



r ...rjv N 



0. 



(57) 



which reflects the N+l elementary transitions with prob- 
abilities Ti within the system, cf. Eq. (|40|) . 



1. Single Level Quantum Dot 

The case N — 1 describes the single resonant level 
model in the large bias limits. For example the waiting 
time distribution wr(z) of a single level quantum dot 
coupled to a left emitter and a right collector reservoir 
with corresponding tunnel rates To — Tl, T± = Tr and 
counting of electrons in the collector only is given by 



wr{z) 



Ti 



z + Trz + Tl 



(58) 



This expression has two poles at z = —Tr and z = — Tl 
and thus, in its pole structure, contains separate infor- 
mation on the two tunnel barriers (left and right). In 



S R (u>) = I 1 - 



2T L T R 
T 2 +uj 2 



I = 



T L T 



L± R 



(59) 



has two poles at w = ±iT which only depend on the 
sum of both tunnel rates and not the individual tunnel 
rates. Thus, already this simple example reveals that 
the physical information contained in the two quantities, 
Sr(u>) and wr(z), is not the same: in order to extract 
both tunnel rates from the noise spectrum, one needs 
additional information such as the absolute value of the 
current /. 

In the time domain, by Laplace back-transforming 
Eq. ([55]) one finds the result first obtained by Davies et 



wr(t) = TrT 



e -T L r _ e -T R r 
' Tr-Tt. : 



(60) 



which has a short-time expansion wr(t) — TrT^t + 
0(r 2 ). The vanishing of wr(t) for r = indicates that 
after tunneling of an electron into the collector the dot is 
in the empty state an no other electron can follow imme- 
diately The corresponding relaxation current Eq. (ffi?)) 



rrclax / 



(t) = I \l - e -( r « +r ^ 4 



(61) 



describes the increase of the current from zero at t = 
(empty dot) towards the stationary current /. 

2. Large Ring N — > oo 

Another interesting case is the ring, Eq. with 
identical rates that are scaled up with increasing N ac- 
cording to 



Ti = (N + 1) 7 . 
The waiting time distributions 



w N {z) = 1 + 



-JY 



(62) 



(63) 



(AT + 1)7 

then yield the same stationary average current I = 1/ (r 



I = I N = l/w' N {Q) =7 



(64) 



for each N, but the Full Counting Statistics following 
from Eq. I0BJ), 



MX) = (JV + 1)7 (e^-l) , 



(65) 



indicates that all the second and higher cumulants van- 
ish for A^ — > 00. The limit of A^ — > 00 in fact leads to 
deterministic transport without fluctuations, i.e. a wait- 
ing time distribution t& oc (z) ee limjv^oo wn{z) = e~ z ^. 
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In the time-domain, the waiting times converge towards 
a delta peak at the inverse transition rate 7, iUoo(t) = 

5 (t — ij , and the 'relaxation' current, i.e. the current 

for empty initial condition at time t = 0, becomes a series 
of delta-peaks, 



1 

n— 

7 



(66) 



cf. Eq. (|26|) . where the term 'relaxation' is of course 
misleading in this limit. The deterministic charac- 
ter of the transport is confirmed by noting that the 
noise spectrum S(u>), Eq. (|2ip . is identical zero since 

i + E± [w 00 (±iu)- 1 -i]~ 1 = Q. 

D. Multi-Level Single Dot 

The next example is a single quantum dot with N levels 
coupled to a left emitter and a right collector via tunnel 
rates 7* (left) and (right), 



/ — 1l Ti ... r w+ i T N 

71 -Ti ... 

72 -r 2 ... 



1L 



V In 
7i + ••• 





■ In- 







(67) 



In the strong Coulomb blockade regime, single electrons 
occupy one of the N levels at a time, and transport is 
from the left to the right with an infinite bias between 
emitter and collector. This is a single-reset system with 
jump operators J a — |a))((a|, a — L/R (left/right), 
where ((R\ = (0, r 1; T 2 , Tjv), \R)) = (1, 0, 0) T , and 
\L)) = 7^ 1 (0, 7l ,7 2 ,...,7 JV ), ({L\ = ( 7L ,0,...,0) T . Note 
that ((R\R)) = ((L\L)) = and therefore J% = j'[ = 0. 
The waiting time distribution for counting only on the 
left or only on the right side, wl{z) — Wr(z) = ib(z), is 
thus 



){ Z ) = {{i\(z-co)- i \i)) = j2^[( z -^r 1 ] 

i=l 

7,r, 



in 



N 



S (z 



1 {z + 1L ){z + v l y 



(68) 



which again follows from considering the ratio of deter- 
minants. In contrast, a splitting C — Cq + Jl + Jr for 
counting on both sides leads to wll(z) — Wrr(z) = 
owing to the single reset character of the system: a tun- 
neling event on the right side can only be followed by a 
tunneling event on the left side, and vice versa. 

The sum over all levels i in Eq. ((68)) reflects the 'paral- 
lel' character of transport through the multi-level system: 
each level contributes with a single waiting time distri- 
bution (analogous to Eq. but with the left tunnel 




FIG. 1: Waiting time distribution w(t) for transport through 
double quantum dot at temperatures T = 4, 8, 16 K at 
e = lmeV, 0.2meV, O.lmeV, -l.OmeV (clockwise). Right 
tunnel rate HTr — 2.5/ieV, other parameters KTl = O.lmeV, 
T c = O.lmeV, electron-phonon coupling parameter g — 0.002, 
cutoff hcuc = 5meV. 



rate Tl replaced by the sum 71, of all left tunnel rates 
7i). For example, if all rates are identical, 7, = T; = 7, 
one is back to an effective single-level case, Eq. with 
= 7 and Fl = Nj which for large N leads to purely 
Poissonian noise properties. 

Some other interesting features of this system can be 
extracted from Eq. (f68|) . The case N = 2 with strongly 
asymmetric rates leads to dynamical channel blockade as 
characterized by huge Fano factors and super-Poissonian 
nois o 18 ' 23 ' 24 : the transport is fast on the time-scale of the 
inverse of the smaller rate, but the electron is occasionally 
trapped in the 'slow' level. The Fano factor F is easily 
obtained via the second moment of the waiting times 2 , 
F = (t 2 )/(t) 2 — 1, and for identical left and right rates 
7; = one has 



F = 



w"(0) 
w'(0) 2 



1 

- 1 = - 
9 



2li 

72 



2 72 
7i 



(69) 



which increases with increasing asymmetry between the 
two rates. 



E. Double Quantum Dot (Strong Coulomb 
Blockade) 

This is the simplest model where quantum coher- 
ence becomes visible in transport^^^i 2 -^! 2 ^!. Spin- 
polarized electrons move between two tunnel-coupled lev- 
els \L) (left) and \R) (right) attached to fermionic reser- 
voirs. The Hamiltonian is a transport version of the spin- 
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boson model (h = 1) 



n 



e 



(70) 



k,a=L,R 

E (^ Q ijO)(a| +J ff. C .) 

k,a=L,R 

\ a -Q + a Q) > U P = ^2 UJ Q a Q a Q' 



Q 



Q 



with pseudo-spin <r z = \L)(L\ - \R){R\, b x = \L)(R\ + 
\R)(L\, the 'empty' state |0), the standard tunnel Hamil- 
tonian for coupling to the reservoirs 7i rcs , and cou- 
pling of the transport electron in the double dot to a 
phonon bath H. p via Ti. ep . One can derive a general- 
ized Master equation in the limit of infinite source-drain 
bia o 25 ' 26 i 27 and in the regime of strong Coulomb block- 
ade, i.e. with only one additional transport electron in 
the double dot. 

The Liouvillian in the basis p = 
Oo, PL, PR, ^Prl,^Prl) has the form 





t 












\ 






T l 











2T C 


c = 















-2T C 









7+ 


-7- 


2 


-7 -£ 




V 





-T c 


T c 


e 


- T -f-l) 



(71) 



with tunnel rates T Q = ^Y,^ \Vg | 2 <5(e-£ fe J, a = L/R, 
assumed as energy-independent rates for electron-phonon 
interactio n 12 ' 28 ! 31 



7 
7± 



A 2 
ttT, 



— + 2T Ae 





-A/a 



4 - -Ae- A /" 



coth 



coth 



£A 
2 



(72) 



A 2 



-A/a 



with a dimensionless coupling constant g, a Dcbyc cutoff 
uj c , the level splitting A = \/e 2 + 4T 2 , and the inverse 
temperature /3 = (fc^T) -1 . These electron-phonon rates 
correspond to a bosonic environment with Ohmic spec- 



tral density p(u>) 



gue 



/l 



'«e(w). 



This is a single-reset system with a jump operator de- 
scribing the tunneling of single electrons from the right 
dot into the collector, J R = \R))((R\ with ((R\ = 
(0, 0, T R , 0, 0) and = (1, 0, 0, 0, 0) T . The relaxation 
current 7 relax (z) and thus, via Eq. (f25|) . the waiting time 
distribution w(z), is given byi 2 - 



■oiax/ \ _ r - / x r_Rr_Lg+(2:) 

1 j - i ^^- z{ [ z + r fi + . 9 _(z)](z + r L ) + (z + r^ + r i ) ff+ (z)} ^ 

g±{z) = 2T C — — 74 
(7 + rfl/2 + z) 2 + 



The stationary current I follows directly from j rolax (,?: — > 
0) ~ //z, or via the stationary state po that is either 
obtained from the eigenvector of C with eigenvalue zero, 
or from solving the linear equation £'p = (1, 0, 0, 0, 0) T , 
where C is obtained by replacing the first row of the 
singular matrix £ with (1, 1, 1, 0, 0), corresponding to the 
normalization condition no + + ur = 1. 

The five poles of w(z) are obtained analytically via 



w{z) 



7 rclax (z) 
1 + / rclax (z) 



(75) 



and in the case of no electron-phonon scattering are ex- 



plicitly given by 



z\ = -Ti 



?2/3 = 
z 4/5 = 

A R = 



r R , .a r 
± * — 1=\ 

2 V2\ 



4T 2 — . 



Tr£ 
Ar 



- 1 



(76) 



One recognizes that the information on the left tunnel 
barrier contained in the pole Z\ = — is completely 
separated from the other poles. This is a consequence of 
the useful relation Eq. (|5ip for an unraveling according 
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to counting electrons both in the left and right leads, 

w RL (z)w LR (z) 



w(z) = W RR (z) + 



1 - w LL 
wrl(z)w lr (z) = w RL (z) 



r, 



z + Tr' 



(77) 



because two subsequent jumps in the emitter (L) or col- 
lector (R) are not possible due to the strong Coulomb 
blockade assumption, i.e. wll = w RR = 0. A jump 
on the right side leaves the system empty, with the triv- 
ial re-charging process as described by wlr{z) = Z+ ^ L 
following, and thus all the relevant information on the 
quantum system is contained in w R l{z) which describes 
the dynamics after a jump into the double dot from the 
left lead. 

In the time domain, it is easier to obtain w(t) directly 
from its definition Eq. ([6]). The result shown in Fig. [T] 
was produced using the matrix exponential in MATH- 
EMATICA. Parameters were chosen close to those used 
in a recent experimental and theoretical analysis of the 
Fano factor in vertical double quantum dots^. The most 
important feature in w(t) is the appearance of oscilla- 
tions with period f» A for A ^> T R , which are due to 
the coherent coupling oc T c between the two quantum 
dots as reflected in the two imaginary parts in the zeroes 
Eq. (|76p . With increasing temperature T of the phonon 
bath, the oscillations become less pronounced although 
they are still visible at relatively large T. 

The short time expansion of w(t) is obtained via 
w(z -> oo) ~ 2T^T R T L z' 4 or directly from Eq. (|2T]) , 



w(t) = ±2T?T r T l t 3 + 0(t 4 ) 



(78) 



which reflects the elementary transitions within the sys- 
tem: tunneling of an electron from the left lead to the 
left dot oc followed by coherent tunneling from the 
left dot to the right dot oc T c 2 , and finally tunneling of 
an electron from the right dot to the right lead oc T^. 
The cubic behaviour, w(r) oc r 3 , for coherent tunneling 
is in contrast to sequential tunneling through N = 2 dots 
(ring example above) where one would find a quadratic 
behaviour of w(t) at small r, cf. Eq. ([57]) . 



F. Example: Three-State System 

As a non-trivial multiple-reset system, consider the 
Anderson single-impurity model in the limit of infinite 
bias, 

Po = -lL-\Pa + li^P] + iRlPl 

P\ = 1l\Pq - (7-RT + r LT)PT + r _RTP2 

Pi = 1LIP0 - {jRl + ?Ll)Pl + r fli p 2 

P2 = -r L TPt - v lipi - (r^T + r fl4 > 2 - (79) 

The system consist of a single electronic level with four 
electronic states: empty (0), spin up/down (|,!) and dou- 
bly occupied (2), and is again coupled to an emitter (left) 



2.5 
2 

1.5 
1 

0.5 



- 






E (three-state) 




r| (three-state) 




E (classical) 




T| (classical) 





0123456789 
Y R /T L (three-state), r L /T R (classical) 
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FIG. 2: Comparison between current entropy E and waiting 
time parameter r\ as a function of "/r/Fl [Fl/Fr] for the 
three-state system Eq. (j 79|) [the classical model Eq. 



and a collector (right). In the infinite voltage limit, rates 
7 atT describe transitions between empty and singly occu- 
pied states with spin a, and rates T aa describe transitions 
between singly and doubly occupied states. Assuming 
spin-independent rates, one can introduce pi = p-^ + py 
and thus in the components p = (po,p\,p2), the Liouvil- 
lian is 



L = 




e iXR j R 
-(r L +7fl) e*x*2T R 
T L e^ -2T R 



(80) 

Xl=Xr=0 



where the counting fields e %XR i L indicate the jump oper- 
ators J R/L . For example, J R = \R1))({R1\ + \R2))((R2\ 
with kcts = (1,0, 0) T , \R2)) = (0,1, 0) T , and 

bras ((Rl\ = (0, 7fl ,0), ((R2\ = (0,21^,0). This is 
thus a multiple (twofold) reset system with J\j h ^ 

and J R i L — 0, which reflects the two possible reset 
states (singly occupied and empty) contributing to trans- 
port. One can check by direct calculation that the FCS 
equation, Eq. P9"| . as obtained from the two by two 
waiting time matrix W(z) (e.g., for counting in the 
collector) coincides with the usual eigenvalue equation, 
det[£(e l *«) - z] = 0. 

The stationary state is obtained as 



Pi 



Pa 



2 — Po, 

1R 



P2 = r^po 

1R TiJ 



1 



7-R 7k T-R 



(81) 



with a stationary current I = 27^,(1 + ^)pa- The (re- 
duced) waiting time distribution for measuring on either 
the left or the right side only, Eq. (|3ip . follows as 



w {r Xz)=wl>{z)=w y n > {z) (82) 
(zT L + 2 lL { lR + T L )) {2Y L T R + -ir{z + 2T r )) 



(z + 2 lL ) ( 7iJ + T L ) (z + j R + T L ) (z + 2T 



R) 
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The short-time limit, Eq. (|40|) , therefore is 



(r) 



(r-0) 



i + r L / 7JJ 



O(r) 



(83) 



which indicates that there is a finite probability to ob- 
serve two electrons just one after another, in contrast to 
single-reset systems. 

On the other hand, the waiting times for counting in 
both leads, corresponding to a splitting C = Co+Jl+Jr, 
lead to 



T) = w LL (0) = w RR (0) 



{1 + Ir/TlY 



(84) 



In contrast to the single-reset systems discussed above, 
this means that there can be two subsequent electron 
jumps on either side. The quantity 77 > offers it- 
self as an experimentally accessible way to quantify the 
multiple-reset character of the system. In fact, rj can be 
compared to the current entropy, Eq. (|53|) . introduced 
above, 



E = Br = Et, = 



\og(\+ lR /T L ) 
1 + 1r/Tl 



+ {L^> R). (85) 



It is interesting to notice that both quantities, 77 and 
E, depend on the same combination of rates and in fact 
display a qualitatively quite similar behaviour, cf. Figj2] 
A further way to quantify the multiple reset character 
of transport is to exploit the fact that the reduced noise 
spectrum, S^ r '(u>), is not simply obtained via the reduced 
waiting time distribution w^ r \z). One can therefore in- 
troduce a 'fidelity' 



S(u) = I 



E[* (r) < 



±iu J y 1 - 1 



(86) 
, (87) 



which is unity for single reset systems and tends towards 
one for multiple-reset system at large frequencies u> as 
a consequence of Eq. (|42|) . The behaviour of F(uj) at 
various right tunnel rates jr = T R for fixed 7l = Tl is 
shown in Fig. [3] The deviation of F(u>) from unity is 
strongest for ^r/Tl « 1. This corresponds well with the 
current entropy maximum in Fig. [2l although the latter 
only depends on the ratio ^r/Vl whereas F(cu) depends 
on all four tunnel rates independently. 

The limit of very large and very small "f R can be 
understood from the waiting time distribution, i.e. 



w(z) 



2Fl 



which corresponds to the empty 



linXyjj^oo ^ v ~, — 21V+Z' 

stationary state with Poissonian waiting time distribu- 
tion, and lim 7fl ^o w{z) = ( 2 rn+z)(r L +z) ' wn i cn corre- 
sponds to a stationary state p — (0, 2T r /T,Tl/T), 
r = Tl + r« of a two-state system (one or two elec- 
trons). In both these limits, F(lo) tends towards unity in 
agreement with Fig. [3] 





1.25 ri 
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1.15 
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1.1 - 




1.05 - 
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0.95 - 
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FIG. 3: Ratio of the spectrum S(u>), Eq. (|86|l . and the proper 
(reduced) noise spectrum S^ r \w), Eq. (|32|l . for the three-state 
system Eq. (|79|l . Frequency uj and tunnel rates are in units 
of Y L = y L . 



G. Classical Transport 

The final example is a system which can accommodate 
an arbitrary large number n of particles, n = 0,1,2... 
with occupation probability p n in the stationary state 
P = (Po>Pi)P2,P3, ■■■)■ Additional single particles enter 
the system at the rate T l from a left reservoir regardless 
of n, and single particles leave the system into the right 
reservoir at rate T^, leading to system transitions n — > 
n — 1 at a rate nT R . The particles are assumed classical, 
i.e. there are no effects due to boson or fermion statistics. 
The Liouvillian is 



L = 




... 
... 

T L 3T R ... 



(88) 



(infinitely many entries). The average number n t of par- 
ticles in the system as a function of time t simply obeys 



dt 



nt 



T R n t . 



(89) 



One decomposes C = Cq + J R + Jl with jump operators 
with matrix elements 

[Jr\u = kV R 5 k ,i-i, [Jl\u = rz,<Sfc,H-f (9°) 

The system has a stationary state po given by the Poisson 
distribution 



a 

Pn = —re 
n! 



(91) 



Using the stationary current I = I R = 1^ = Ti and the 
stationary state po with Eq. (|9ip . the waiting time dis- 
tributions are obtained by simple matrix multiplication 
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in uiij(z) = Tr JiWj(z)po/I, 
wll(z) = w RR (z) 

wlr(z) 
wrl(z) 



apr, 



^— ' z/T R + a + n + 1 

n— 



t'o z I Tr + 01 + n - 1 



(n + l)p v 



oo 

^ z/T R + a + n + 1 

n— 



, (92) 



which can be expressed in terms of (incomplete) Gamma 
functions and which correspond to a measurement on 
both the right and the left side. A measurement on, e.g., 
the left side only corresponds to a splitting C = £q + Jl 
with a single jump operator J7l only. The corresponding 
waiting time distribution, wl{z), has to fulfill the iden- 
tity Eq. (fSTj) between the two splittings from Sect. IIII A[ 



w L (z) = w LL (z) + 



w LR (z)w RL (z) 



1 - w RR (z) 

which by explicit calculation using Eq. (f9"2")) yields 



w L (z) 



z + r L ' 



(93) 



(94) 



which indeed is the expected result as the particles on the 
left side enter the system independently. Since wll(z) = 
w RR {z), one obtains the same result, wl(z) = w R (z), for 
counting on the right side. 

There are an infinite number of reset states with cor- 
responding currents IkL = Pk^L for k = 0,1,2,... and 
hit = PkkT R for k = 1,2,.... Using kp k = apk-i, 
Eq. (|9"Tj) , the current entropies, Eq. ([55)1 . are given by 
the entropy of the Poisson distribution p n itself, 



E = E R = E L = S[ Pn ] = - J^Pn lo. 



gPr, 



(95) 



which is positive for a = Tl/T r > 0. On the other hand, 
the quantity 



r) = wll(0) = w RR (0) 



(96) 



is non-zero as there is a finite probability for yet another 
transition into or out of the system immediately follow- 
ing a previous one. Both the entropy E and the zero- 
frequency waiting time r\ thus reflect the multiple-reset 



character of the system. As a function of the asymmetry 
parameter a, both are functions that increase from at 
a = with growing a, cf. Fig. [21 but 77 has an asymp- 
tote at \ whereas E continues to grow logarithmically 32 
as E ~ i log(27reo;). 



IV. REMARKS AND CONCLUSION 

Other multiple-reset situations include systems with 
internal degrees of freedom. Transport of electrons 
through a single level with strong coupling to a vibra- 
tional degrees of freedom leads to a block structure of 
the Liouvillian: the simple (scalar) entries in the two-by- 
two matrix of the N = 1 case, Eq. ([54]) . become matrices, 
where the relatively complicated structure of the matrix 
elements leads to an avalanche- type of transport with 
non-trivial power-laws in the noise spectrum S(ui)^. It 
would be interesting to find other systems where similar 
features can be extracted from waiting times of Liouvil- 
lians with a relatively complex structure. 

As a conclusion, the waiting times w(t) appear to be 
a flexible theoretical tool for describing single particle 
transport, in particular as they contain the other statis- 
tical quantities - FCS and noise spectrum S(u>) - that 
have mainly been used in the past. Depending on the 
system, measuring w(t) could provide additional infor- 
mation or a least serve as a cross-check for FCS and noise 
data. 

One interesting generalization of the formalism in Sect. 
II would be to consider Master equations of integro- 
differential type, e.g. containing non-Markovian mem- 
ory kernels and retardation effects due to, e.g., strong 
electron-phonon couplin g 33 i 34 . A further question is the 
connection between un-symmetrized noise spectral and 
w(t), and the relevance of multi-time waiting time distri- 
butions (n > 2 jump operators in Eq. Q) and their re- 
lation to higher order frequency-dependent cumulants 22 . 
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